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ABSTRACT 

Results from the first dedicated study of Galactic planetary nebulae (PNe) by means of op- 
tical integral field spectroscopy with the VLT FLAMES Argus integral field unit (IFU) are 
presented. Three typical Galactic-disk PNe have been mapped with the 11. "5 x 7. "2 Ar- 
gus array: two dimensional spectral maps of the main shell of NGC5882 and of large areas 
of NGC6153 and 7009 with 297 spatial pixels per target were obtained at sub-arcsec res- 
olutions. A corresponding number of 297 spectra per target were obtained in the 396.4 - 
507.8 nm range. Spatially resolved maps of emission lines and of nebular physical properties 
such as electron temperatures, densities and ionic abundances were produced. The abundances 
of helium and of doubly ionized carbon and oxygen, relative to hydrogen, were derived from 
optical recombination lines (ORLs), while those of 2+ were also derived from the classic 
collisionally excited lines (CELs). The occurrence of the abundance discrepancy problem, 
pertaining to oxygen, was investigated by mapping the ratio of ORL/CEL abundances for 
2+ (the abundance discrepancy factor, ADF) across the face of the PNe. The ADF varies 
between targets and also with position within the targets, attaining values of ~ 40 in the case 
of NGC6153 and ~ 30 in the case of NGC7009. Correlations of the ADF with geometric 
distance from the central star and plasma surface brightness (for NGC 6153), as well as with 
[O III] electron temperature, plasma ionization state and other physical properties of the tar- 
gets are established. Very small values of the temperature fluctuation parameter in the plane 
of the sky, t\(0 2+ ), are found in all cases. 

It is argued that these results provide further evidence for the existence in run-of-the- 
mill PNe of a distinct nebular component consisting of hydrogen-deficient, super-metal-rich 
plasma. The zones containing this posited component appear as undulations in the C II and 
O II ORL abundance diagnostics of about 2 spatial pixels across, and so any associated struc- 
tures should have physical sizes of less than ~ 1000 astronomical units. Regarding the origin 
of the inferred zones, we propose that circumstellar disks, Abell 30-type knots, or Helix-type 
cometary globules may be involved. Implications for emission line studies of nebulae are 
discussed. 
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1 INTRODUCTION 

Planetary nebulae represent only an ephemeral stage in the late 
evolution of low- to intermediate-mass stars, the predominant stel- 
lar population in late type galaxies, but their rich emission line 
spectra allow us to quantify the amounts of helium, and heav- 
ier elements (mainly carbon, nitrogen, and s-process elements) 
that these stars produce and return into the interstellar medium 
(ISM). Recent studies have indicated that PNe may also contain 
non-negligible amounts of endogenous oxygen and neon synthe- 



sized by intrinsic nucleosynthesis in the stellar progenitors and ex- 
pelled into the nebulae following the third dredge-up, especially 
more so in lower metallicity environments such as those of the 
Galactic Halo (Pequignot & Tsamis 2005) and of satellites of the 
Galaxy (Pequignot et al. 2000; Zijlstra et al. 2006; Leisy & Den- 
nefeld 2006). Their importance as vital sources of neutron-capture 
elements, such as for instance selenium and krypton among oth- 
ers, is now observationally confirmed (Pequignot & Baluteau 1994; 
Sterling et al. 2007). With the advent of 8-m class telescopes, PNe 
are increasingly being used as probes of the chemical evolution of 



2 Y. G. Tsamis et al. 



Local Group galaxies (e.g., Goncalves et al. 2007; Saviane et al. 
2008), while by virtue of their intrinsically bright monochromatic 
[O III] 500.7 nm emission they have been detected in the intraclus- 
ter stellar population of the Coma cluster at a distance of 100 Mpc 
(Gerhard et al. 2005, 2007). 

A focal problem, however, remains in the astrophysics of these 
objects: namely that large discrepancies are measured between 
abundances and temperatures obtained from recombination lines 
and continua versus those obtained from the bright forbidden lines 
(e.g., Tsamis et al. 2003b, 2004; Liu et al. 2000, 2004; Wesson et 
al. 2005). For the majority of more than 100 PNe surveyed thus far 
by long-slit (or echelle) spectroscopy, abundances of the elements 
carbon, oxygen, nitrogen and neon derived from their optical re- 
combination lines (ORLs) are higher than those derived from the 
classic ('forbidden') collisionally excited lines (CELs) by factors 
2-3 (see review by Liu 2006 and references therein). For about 5- 
10 per cent of the surveyed objects however the discrepancies are 
in the 4-80 range, with the most pathological case to date being 
Hf2-2 (Liu et al. 2006), a southern planetary nebula. In the cases 
of carbon, nitrogen (Tsamis 2002; Tsamis et al. 2004) and oxygen 
(Liu et al 2001; Tsamis et al. 2004), large metal overabundances are 
correlated with lower plasma temperatures derived from ORLs and 
continua than from forbidden lines. The resolution of this problem 
remains of high priority as it continues to undermine our trust on 
secure elemental abundance determinations for PNe, and for other 
types of emission line objects (such as H II regions; e.g. Tsamis et 
al 2003a). 

In the overwhelming majority of cases the ratio of ORL/CEL 
abundances for a given heavy element ion (the abundance discrep- 
ancy factor, ADF) is not correlated with the excitation energies, 
Sex, of the CELs involved: meaning that, per individual PN, in- 
frared (IR) or ultraviolet (UV)/optical CELs of low and high E czi 
respectively, yield very comparable ionic abundances, which are 
uniformly lower than the corresponding abundances of the same 
ions derived from ORLs. The most likely explanation for these 
spectroscopic results was judged to be the presence of cold plasma 
regions embedded in the nebular gas in the form of relatively 
dense, hydrogen-poor condensations - clumps or filaments (Liu 
et al. 2000; Tsamis 2002, Tsamis et al. 2004). Due to their ele- 
vated content in heavy elements (several times Solar), these plasma 
regions would have reached photoionization equilibrium at lower 
temperatures than the ambient 'normal' composition gas by emit- 
ting far-IR lines, the primary nebular thermostat in relatively low 
temperatures. Since the emissivity of heavy element recombination 
lines is enhanced at lower temperatures, while at the same time 
that of the classic forbidden lines is diminished, the hydrogen-poor, 
heavy element-rich clumps would emit heavy element recombina- 
tion lines profusely, yielding a truer estimate of the heavy element 
content of these regions. 

Whereas the long-slit surveys yielded concrete evidence for 
elevated recombination-line abundances, the temperature of the 
suspected hydrogen-deficient clumps was tougher to determine 
since the diagnostic lines involved are faint and often suffer from 
blends in lower resolution spectra. Nevertheless, there has been ev- 
idence from those studies for temperatures lower by several thou- 
sand K than the typical CEL temperatures of photoionized neb- 
ulae (Tsamis et al 2004; Liu et al. 2004). This is one aspect of 
the proposed 'dual abundance model' solution that renders it self- 
consistent, something rather lacking from alternative propositions, 
such as small-scale temperature fluctuations in a chemically homo- 
geneous medium (e.g. Peimbert et al. 2004), which do not invoke a 
physical process to explain the discrepancies. A major unresolved 



issue remains the as yet unknown origin of the high abundance 
clumps, and their mass relative to the total ionized nebular gas. 

Heavy element recombination lines in nebulae can be a thou- 
sand times or more fainter than hydrogen recombination lines, in 
sharp contrast to the luminous forbidden lines emitted by the same 
heavy ions. We therefore brought to bear on this thorny issue the 
combination of integral field spectroscopy (IFS) and the light col- 
lecting power of the Very Large Telescope (VLT). IFS techniques 
allow spectra to be obtained from every spatial element in a two- 
dimensional field of view, combined with imaging capability of 
the field at any wavelength across a given A range. This is usually 
achieved by means of an IFU array composed of microlenses cou- 
pled to optical fibres or slicers (e.g. Allington-Smith 2006). These 
are fed into a spectrograph and are subsequently individually im- 
aged on the focal plane of the telescope. 

The contiguous array of microlenses (or spatial pixels - 'spax- 
els') can also be thought of as a series of long-slits stacked together 
and, in this sense, an IFU ofm x n spaxels used on a single tele- 
scope is equivalent to using an array of m telescopes equipped with 
long slits of n pixels each. It could be argued that in the former case 
a given observing programme is executed m times faster than when 
using a classic long-slit spectrograph. In reality this is not strictly 
true as invariably there are some light losses associated with IFU fi- 
bres and their optical connection to the spectrograph (particularly at 
the interfaces with the microlenses and the feed into the slit). Over- 
all, IFU observations afford considerable advantages especially in 
a science programme like ours where faint diagnostic lines are tar- 
geted, and where 2D mapping of nebular regions is sought in the 
hope of isolating suspect metal-rich regions which may a priori 
have a random spatial distribution in the nebula. 

In addition, IFS data can provide the much-needed input to 
fully realize the capabilities of new plasma codes such as those 
developed in recent years for the 3D modelling of photoionized 
nebulae (e.g. Ercolano et al. 2003; Morisset et al. 2005). For ex- 
ample, such codes can produce simulated 3D or 2D spectral im- 
ages of model nebulae at arbitrary wavelengths and inclinations 
that can be readily evaluated using images obtained by IFS. More- 
over, high spectral resolution IFS data such as those we obtained 
for NGC 7009 (Walsh et al. in preparation) can yield detailed in- 
formation on the velocity field of an object that could allow a 3D 
kinematical model to be built. These are some of the wider aspects 
of this work. 

Here we pursue a physicochemical analysis of three repre- 
sentative PNe belonging to the Galactic-disk population with pub- 
lished oxygen ADFs ranging from ~ 2 - 10. Long-slit spectroscopy 
of the targets has been published: NGC 5882 (Tsamis et al. 2003b, 
2004); NGC 6153 (Liu et al. 2000); NGC 7009 (Liu et al. 1995; 
hereafter LSBC). In those studies, the first two targets were ob- 
served with a long slit scanned across the whole nebular surface. 
This technique, which was applied throughout the long-slit surveys, 
yielded mean relative emission line fluxes for the whole PNe, but 
naturally any variations of the physical properties in the direction 
perpendicular to the slit could not be sampled. An analysis along a 
NGC 6153 fixed slit sampling a 'slice' of the nebula revealed sub- 
stantial variations of the physical properties over the slit length. 
Typically, however, the formerly employed ~ 2 arcsec-wide slits 
achieved a rather coarse spatial sampling of ~ 1.6 arcsec per pixel, 
without sampling the seeing which was of the order of 1 arcsec. 

We obtained contiguous 2D spatial coverage of these PNe at 
sub-arcsec resolution with the Fibre Large Array Multi Element 
Spectrograph (FLAMES )/Giraffe IFU Argus. Our main aim was 
to map out the nebular properties and their variations across large 
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Figure 1. NGC 5882: VLT Argus H/3 map overlaid on a HST WFPC2 image 
of the nebula covered by the 11. "5 X 7. "2 IFU. The HST snapshot was 
taken with a broad-band F555W filter (credit: H. Bond and Space Telescope 
Science Institute). The position of the stellar nucleus is visible in the inset. 




Figure 2. NGC 6153: VLT Argus H/3 map overlaid on a pseudocolour HST 
WFPC2 broad-band filter F814W image of the nebula covered by the 11. "5 
X 7. "2 IFU (source of Hubble image: Liu et al. 2000). The position of the 
stellar nucleus is visible in the inset. 



nebular areas and to investigate in detail the existence of any heavy 
element-rich nebular regions, so as to place constraints on models 
that have been put forward to explain this astrophysical quandary. 
The sample PNe were chosen because NGC 6153 and NGC 7009 
exhibit the abundance anomaly problem to a high degree (with 
known oxygen ADFs of ~ 10 and 5 respectively) and are excellent 
cases for a detailed study, whereas NGC 5882 (with an ADF of ~ 
2) was meant to be our 'control' target as it represents the majority 
of Galactic nebulae that show a moderate degree of abundance dis- 
crepancies. The paper layout is as follows: the specifications of the 
observational data set are discussed in Section 2, the data reduction 
methods are presented in Section 3, and the results in Sections 4 
and 5. The discussion and our conclusions appear in Section 6. An 
introduction to this study was presented by Tsamis et al. (2007). 



2 OBSERVATIONS 

Integral field spectra of NGC 5882, 6153 and 7009 were obtained in 
service mode between March-June 2005 in sub-arcsec seeing with 
the FLAMES/Giraffe Argus IFU on the 8.2-m VLT UT2/Kueyen 
(P.I. Tsamis; ESO Period 75). The observing log is presented in Ta- 
ble 1. All targets were observed in the 396.4 - 507.8 nm range with 
the large 11. "5 x 7. "2 IFU. The Argus unit is a rectangular array of 
22 x 14 microlenses fed by optical fibres: the 11. "5 x 7. "2 array 
thus projects 308 spaxels on the sky; the individual spaxel size is 
0."52 x 0."52. Eight spaxels are reserved for sky subtraction pur- 
poses (from a total of 14 dedicated sky-fibres) and three currently 
correspond to broken fibres (appearing in the second row of a re- 
constructed Argus image; see Fig. 6); the remainder provide 297 
individual spectra per target. The Giraffe spectrograph is equipped 
with a 2048 x 4096 EEV CCD which has 15 /zm pixels. Two sep- 
arate wavelength ranges were observed using the LR02 (3964— 
4567 A) and LR03 (4501-5078 A) grating set-ups. The respective 
spectral resolving powers, R (= A/AA), were 10 200 and 12 000. 
Higher spectral resolution HR04 and HR06 spectra (R = 32 500) 
of NGC 7009 were also taken with the small 6. "6 x 4."2 IFU 
and will be published in a separate paper (Walsh et al. in prepara- 
tion). Total exposure times were 6930 sec (NGC 5882), 24 000 sec 
(NGC 6153), and 4260 sec (NGC 7009); more time was devoted to 
NGC 6153 than to the other two targets since this nebula has the 
lowest surface brightness, and exhibits the abundance anomaly to 
a high degree. The average seeing during the NGC 5882 exposures 
was 0.6 arcsec full width at half maximum (FWHM); during six 
out of ten observing blocks of NGC 6153 it was better than 0.8 arc- 
sec, while during the remaining four it was ~ 0.9 arcsec. During 
the NGC 7009 exposures the seeing ranged between 0.5-0.8 arc- 
sec. The sky transparency conditions were generally in the 'clear' 
category. 

The positions observed on each PN were given due consid- 
eration. In Figs. 1-3 we show the exact positions observed by 
overlaying Argus spectral images taken in the light of H/3 4861 A 
on Hubble Space Telescope (HST) WFPC2 images of the nebulae. 
The Argus maps have been drawn to scale in all cases. Regarding 
NGC 5882 (Fig. 1), which is the nebula with the smallest apparent 
size of the three, the IFU array was positioned so that the largest 
part of the bright shell of the PN as well as part of the fainter halo 
to the south could be covered. The centre of the IFU was positioned 
at a = 15 h 16 m 49 s . 8, 5 = -45° 39' 01".4 (J2000). 

For NGC 6153 (Fig. 2) we chose to align the large axis of 
Argus roughly along the minor axis of the PN so as to cover both the 
bright patch visible on the HST image in the southeastern outskirts 
of the nebula and the inner fainter regions closer to the central star 
- this position also shares the orientation of the 2-arcsec wide fixed 
long slit employed by Liu et al. (2000). The centre of the IFU was 
positioned at a = 16 h 31 m 31 s .O, 5 = -40° 15' 12" .4 (J2000). 

NGC 7009 (Fig. 3) has an elliptical shape with low ionization 
emission regions at both ends of its major axis of symmetry which 
give it a bipolar appearance. It shows a clearly defined bright in- 
ner shell surrounded by a fainter outer envelope. We thus chose to 
place the IFU array so that both the faint regions in the centre of 
the PN and the sharp rim of the bright inner shell could be cov- 
ered. The centre of the IFU was positioned at a = 21 h 04 m 10 s . 8, 
S = -11°21'48".7 (J2000). In all cases the central stars of the 
PNe (visible in the insets of the HST/Aigus overlays) fell within 
the IFU field of view and were a useful point of reference. This fur- 
ther allowed us to examine the nebular physical conditions in the 
immediate vicinity of the stellar nuclei. 
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Table 1. Journal of VLT 8.2-m FLAMES Argus observations. 



Target 


Date 


A-range 


Grating 


Exp. time 




(UT) 


(A) 




(sec) 


NGC5882 


23/03/05 


3964^1567 


LR02 


3x440 




23/03/05 


4501-5078 


LR03 


3x330 




29/04/05 


3964^4567 


LR02 


6x440 




29/04/05 


4501-5078 


LR03 


6x330 


NGC6153 


29/04/05 


3964-4567 


LR02 


6x800 




29/04/05 


4501-5078 


LR03 


6x400 




30/04/05 


3964^4567 


LR02 


12x800 




30/04/05 


4501-5078 


LR03 


12x400 




01/05/05 


3964-4567 


LR02 


4x800 




01/05/05 


4501-5078 


LR03 


4x400 


NGC7009 


02/06/05 


3964-4567 


LR02 


8x270 




02/06/05 


4501-5078 


LR03 


12x175 




Figure 3. NGC 7009: VLT Argus H/3 map overlaid on a HST WFPC2 im- 
age outlining the area of the nebula covered by the 11. "5 X 7. "2 IFU. The 
//Srpseudocolour is a composite of F658N = [N II] (red), F502N = [O III] 
(green) and F469N = He II (blue) frames (credit: B. Balick and Space Tele- 
scope Science Institute). The position of the stellar nucleus is seen in the 
inset. 



3 DATA REDUCTION 

The spectra were reduced at ESO Garching with the girBLDRS 
pipelinf] in a process which included bias removal, localization 
of fibres on the flats, extraction of individual fibres, wavelength 
calibration and rebinning of the ThAr lamp exposures, and the 
full processing of the science frames which resulted in cosmic- 
ray cleaned, flat-fielded, wavelength-rebinned spectra. In practice 
the pipeline actions were initiated by the creation of a master bias 
frame from the raw biases (BIASMAST pipe) which was subtracted 
from the science exposures at a later reduction step. The localiza- 
tion of fibres on the flat field frames followed (LOCMAST pipe). 
The full master wavelength solution was produced via the extrac- 
tion, wavelength calibration and rebinning of the raw Th-Ar expo- 
sures (WCALMAST pipe). This was followed by the full processing 
of the science frames via the EXTRACT recipe which, by using the 
products of the previous steps, subtracts the master bias from the 
science exposures, takes out an average rejecting any cosmic ray 



The girBLDRS pipeline is provided by the Geneva Observatory 
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Figure 4. Spectrum of NGC 5882 from spaxel (11, 13) highlighting the O II 
VI multiplet recombination lines. The red line is a sum of Gaussians plus a 
cubic spline continuum fit to the data (black line). The top panel shows the 
residuals of the fit (black) compared to the error bars of each point in the 
spectrum (red). 



hits, performs the flat-fielding, and finally rebins the extracted spec- 
tra in wavelength space. In this way, four to six science frames per 
observing block were processed together and the resulting frame 
was corrected for atmospheric extinction within IRAF. The airmass 
during the exposures was ~ 1.1. When this varied by more than 
~ 0. 1 during the course of a single observing block (1 hr) the 
frames were extracted individually and were corrected for atmo- 
spheric extinction one at a time. They were subsequently averaged 
using the IMCOMBINE task which performed the cosmic ray re- 
jection as well. The flux calibration was achieved using exposures 
of the white dwarf EG 274 using the tasks CALIBRATE and STAN- 
DARD. The flux standard exposures were first individually extracted 
with girBLDRS and the spaxels containing the star were summed 
up to form a ID spectrum. The sensitivity function was determined 
using IRAF's SENSFUNC and this was subsequently applied to the 
combined science exposures. The sky subtraction was performed 
by averaging the spectra recorded by the sky fibres and subtract- 
ing this spectrum from that of each spaxel in the IFU. Each of the 
14 sky-fibre spectra were examined and the average count in each 
compared. Four sky fibres had higher counts caused by contamina- 
tion from the simultaneous calibration lamp (see below), and were 
not included in the computation of the mean sky spectrum. 

Custom-made scripts allowed us to convert the row by row 
stacked, processed CCD spectra to data cubes. The nebular emis- 
sion line fluxes were generally measured automatically by fitting 
Gaussian profiles via means of the dedicated program SPECTRE 
which is based on MINUIT \ 2 minimization routines (James 1998). 
This enabled us to retrieve from the spectra associated with each 
spaxel: the line fluxes, the line FWHM, equivalent widths (EW), 
line centroids, and the respective errors of these quantities which 
were propagated in all physical properties derived afterwards. 
The same program was used to interactively fit the nebular lines 
in the vicinity of the PN nuclei, and the absorption lines under- 
neath the H I lines over the same spaxels. In Figs. 4 and 5 Giraffe 



2 The pre-analysis errors produced by the pipeline were used as the starting 
point of the error analysis; they also include uncertainties from the flat- 
fielding and the geometric correction of the spectra and not just photon- 
noise (Blecha & Simond 2004). 
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Figure 5. Same as Fig. 4 for NGC 6153 (top) and NGC 7009 (bottom). 



spectra of the PNe from single spaxel extractions in the vicinity 
of 465.0 nm are shown, covering the full O II VI recombination 
multiplet, a few N II and C III ORLs, N III Bowen fluorescence 
lines, and He II A4686. The lower panels display the observed spec- 
trum (black line), along with the sum of the fitted continuum (cubic 
spline fit) and the individual Gaussians (red line). The upper pan- 
els show the errors bars at each point in the spectrum (red) and the 
residuals from the line and continuum fitting (black). In general the 
fit residuals are comparable to the error bars. The fits in Figures 
4-5 are for single Gaussian line fits to each emission species; for 
NGC 6153, 130 spaxels were fitted with double Gaussians, since 
the velocity splitting (due to nebular expansion) of 0.6 A could be 
resolved by two lines of FWHM 0.5 A. In the case of NGC 5882, 
30 spaxels were fitted by double profiles. 

The Gaussian fits to the brightest lines in Figures 4-5 show 
that the detailed profile deviates from a Gaussian and in the case 
of NGC 7009, for example, the very strong He II A4686 line shows 
residuals in the line wings that are greater than the errors. Since 
only Gaussians were used for fitting, without a detailed line pro- 
file describing the grating scattering line wings, these line wing 
features were not fitted in detail. Neglect of these faint features 
constitutes a loss of less than one per cent to the line flux, and it 
was considered that a consistent line profile shape for all lines (i.e. 
Gaussian) was of highest priority for comparison of relative line 
fluxes (on which plasma conditions and abundance determinations 
that are central to this analysis depend). In general, the spectral res- 
olution achieved was sufficient to resolve any known blends that 
affect heavy element ORLs in lower resolution spectra. Extra care 
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Figure 6. A spectral map of the continuum showing the characteristic pat- 
tern caused by contamination of the CCD from light leaking from the si- 
multaneous calibration lamp of FLAMES. See text for details. 



was taken to interactively fit any remaining blended lines and to 
also fit double Gaussian profiles in cases where the spectral resolu- 
tion was high enough to resolve the expansion of the nebula (e.g. 
in the case of NGC 6153). 

An unforeseen issue arose due to the use of the 'simultaneous 
calibration lamp' during the science exposures. Use of this device is 
not really necessary for spectrophotometric work and it was left on 
unintentionally during the preparation of phase II. The lamp intro- 
duced an additive contamination pattern (Francesca Primas, private 
communication) on certain spaxels of the array (see Fig. 6), which 
showed up as a mix of broad and narrow features. All these features 
are light scattered from fibre-to-fibre and the broad features are 
presumably the scattering wings of strong lines in the comparison 
lamp spectrum. Since the PN line emission covered the whole of 
the IFU area for all three nebulae, there was no opportunity with the 
current data to subtract PN emission-free spaxel spectra in order to 
remove the contamination pattern. Instead the contaminated spectra 
were interactively fitted by following the continuum in detail and 
any spurious narrower 'lines' were fitted by additional Gaussians. 
The results were satisfactory, although of course the fit errors on 
the nebular line fluxes were increased at those spaxels. The spectra 
of NGC 6153, which is the nebula with the lowest surface bright- 
ness, were most affected and 40 spaxels were interactively fitted 
and decontaminated in the manner described above. 

Likewise, a scattered light feature appeared towards the last 
few CCD columns (long wavelengths) at the top-right corner of 
the stacked array. This was due to a residual CCD signal, which 
is a known cosmetic defect of the chip. The signal is additive (in 
excess of the bias) and, as it affected only certain spaxels in the re- 
constructed Argus image, it was fitted with a low-order polynomial 
and subtracted out. Finally, apart from the three fibres which were 
broken during ESO Period 75, the fibre corresponding to spaxel 
(4, 4jf| of the reconstructed array is problematic in that it exhibits 
an anomalous spike in its signal, which often affected spaxels (3, 
4) and (5, 4) too - the respective spectra were therefore handled 
accordingly and generally suffer from large uncertainties. 



(X, Y) spaxel coordinates in the text are measured respectively on the 
minor and major axis of the Argus IFU; see Fig 6. 
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4 RESULTS 

After the line fluxes had been measured, the reconstructed Argus 
data cubes were corrected for interstellar extinction using the c(H/3) 
extinction maps derived from a comparison of the observed to the 
predicted H7/H/3 line ratio using its case B value from Storey & 
Hummer (1995). The Galactic extinction law of Howarth (1983) 
with a ratio of total to selective extinction (Rv) of 3.1 was used. 
Since the LR02 and LR03 spectra, containing the H7 and H/3 lines 
respectively, were taken in different exposures, and often on differ- 
ent nights, then it was necessary to tie the fluxes together in order to 
determine the extinction of the maps from the H7 to H/3 ratio. This 
comparison was performed by taking the total flux from the H7 and 
H/3 maps and comparing the ratio with that observed from long-slit 
studies. For NGC5882 the absolute value of the ratio was taken 
as 0.410 from the scanned-slit (mean over the whole PN) data of 
Tsamis et al. (2003b); for NGC6153, the ratio for the whole neb- 
ula (0.333) was taken from Liu et al. (2000); for NGC7009, the 
ratio of 0.443 for a value of c(H/3) = 0.20 was adopted (cf. LSBC). 
It should be noted that the IFU area does not exactly correspond 
to the comparison regions from the long-slit data, so the absolute 
extinction, and hence dereddened line ratios can be systematically 
offset. To gain an idea of the photometric quality of the observa- 
tions, the LR02 and LR03 maps had to be corrected by 47, 5 and 5 
per cent for NGC 5882, 6153 and 7009 respectively. 

Potential shifts between the LR02 and LR03 maps, caused ei- 
ther by a slightly different pointing between the observations (al- 
though note that the LR02 and LR03 setting were always taken in 
pairs; Table 1) and/or the changing differential atmospheric extinc- 
tion between the settings, were determined by cross-correlation of 
features in the H7 and H/3 maps. For NGC 5882 a shift of 0.25 
spaxels in the IFU X-direction could be clearly detected (this PN 
has the sharpest spatial features as it was observed under ~ 0.6 
arcsec seeing). It was nonetheless decided not to correct the maps 
for this small shift since data at the IFU edges would be compro- 
mised. The extinction map for NGC 5882 does not show a very 
pronounced cusp at the position of the shell, as might be expected 
if the H7 and H/3 maps were poorly aligned. 

In the paragraphs that follow we present spectral maps for rep- 
resentative emission lines, the plasma physical conditions which 
were based on spectral map ratios, and the resulting chemical abun- 
dances. 

4.1 Emission line maps 

4.1.1 NGC 5882 

In Fig. 7 a set of NGC 5882 maps is presented showing the ob- 
served H/3, the logarithmic extinction at H/3, and a sample of rep- 
resentative dereddened spectral maps in the light of He I A4471, 
He II A4686, [O III] A4959, [Ne III] A3967, and the recombination 
lines C II A4267 and O II A4649. The top two rows of the array ap- 
pear very bright in the 'unmasked' H/3 image: this does not appear 
in images of the other two PNe [except perhaps in spaxels (3-5, 22) 
for the brightest emission lines of NGC6153] and was judged to be 
a combination of the effects of differential atmospheric extinction 
at the edge of the field of view and the fact that the line flux there 
shows quite a steep gradient - this is also evident from a compari- 
son of our line map with the HST F555W image of the PN which 
shows a confluence of clumped, relatively high surface brightness 
emission in that area (Corradi et al. 2000). Such a gradient at the 
edge of the array could introduce errors in ratios of lines originating 
from the different LR02 and LR03 gratings and so these two rows 



were masked out and are not considered in the physical analysis of 
NGC 5882 that follows. We also decided to exclude the faint PN 
halo from a detailed analysis as the line S/N ratio there, especially 
for the O II ORLs, was much lower than the IFU mean. A mask 
based on the H/3 signal-to-noise ratio was therefore applied to ex- 
clude the halo region which appears with zero value spaxel intensi- 
ties (coloured deep blue) in the F(H/3) map of Fig. 7. Mean values 
of the physical conditions in the halo are given in Sections 4.2. 1 
and 4.3.1. 

The emission lines peak on the main shell of NGC 5882 and 
especially on the north-west corner which coincides with the bright 
patch visible on the HST image (see inset of Fig. 1), our standard 
measure of comparison. The PN appears elliptical as in the Hubble 
image, and the southern portion of its halo is clearly registered in 
the 'unmasked' H/3 map. The south segment of its main shell ap- 
pears faint in all lines and coincides with a low surface brightness 
blister-like feature on the HST image which extends towards the 
halo. The central regions are also much fainter relative to the main 
shell. The [O III] and [Ne III] maps appear very similar, as do the 
C II and O II ORL maps. He II A4686 is more concentrated on the 
west part of the PN than on the east side, in contrast to the He I 
A4471 emission which traces almost the full nebular ring. The dust 
extinction is variable across the nebula; c(H/3) has a mean value of 
0.449 ± 0.034 with an rms deviation of 0.241. Assuming a uniform 
interstellar extinction screen over the face of the PN, the central- 
eastern regions appear to be less dusty than the nebular shell. 



4.1.2 NGC 6153 

In Fig. 8 a set of NGC 6153 maps is presented showing the ob- 
served H/3, the logarithmic extinction at H/3, and a sample of rep- 
resentative dereddened spectral maps of metallic CELs, the helium 
recombination lines and the C II and O II ORLs. The Argus IFU tar- 
geted the brightest portion of the south-east quadrant of the nebula. 
The most striking feature in the maps is the high surface bright- 
ness patch peaking at around spaxel (8, 17) and which is clearly 
identified with the bright area in the south-eastern quadrant of the 
HST image (cf. Fig. 2). This area seems to have lower than average 
dust extinction associated with it, at c(H/3) ~ 0.9 - the extinction 
is relatively higher towards the inner part of the nebula, and in the 
vicinity of the central star, with values closer to 1.3. The variations 
may not be very significant statistically; c(H/3) has a mean value 
of 1.146 ±0.180 with an rms deviation of 0.120. A lower surface 
brightness discontinuity is present in the bright area at coordinates 
~ (5, 14) in all emission lines except [Ar IV] A4711 (which appears 
bright over a large region). This discontinuity can be identified with 
the 'thinning' of the bright patch towards the north-east, something 
also seen in the inset of Fig. 2. He I A4471 covers the same area of 
the IFU as H/3. The He II A4686 peak does not coincide with the 
H/3-bright region but lies interior to it, by 2 spaxels, towards the 
inner nebula where the plasma ionization state is higher - there is 
also a local He II enhancement near the PN nucleus. The [Ne III] 
and [O III] maps are very similar and the lines peak at the same 
spots. The C II and O II ORL maps are also very similar and show 
peaks near the H/3 maximum, but also local peaks in the vicinity of 
the central star. Spaxels (3-5, 22) for the brightest emission lines of 
this PN appear with higher values than the mean, and may suffer to 
some unspecified degree from the same problem as row Y = 22 of 
NGC 5882. At the same time, the inset of Fig. 2 shows that there is 
a local flux enhancement within that IFU area. 
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Figure 7. NGC5882 monochromatic line flux (F) and dereddened intensity (I) maps (in units of 10 — 17 erg cm" 2 s" 1 ), along with the extinction constant 
employed, c(H/3). Blank spaxels: three in 2nd row from top (broken Argus fibres); two per corner (fibres reserved for sky subtraction). The top two rows have 
not been considered in the physical analysis - see text for details. The central star position is at coordinates (6, 15). North is up and east is to the right (cf. 
Fig. 1). 



4.1.3 NGC7009 

In Fig. 9 we present a similar set of maps for NGC7009 as for 
NGC6153. The Argus array targeted the northern part of the in- 
ner shell of the nebula and the long IFU axis lies in approximately 
the north-east to south-west direction (cf. Fig. 3). The most con- 
spicuous feature in all maps is the bright arc that runs from top 
to bottom and shows high surface brightness spots at either end - 
this is mostly evident in the light of H/3 but is also present in other 



wavelengths too. It can be safely identified with the rim of the inner 
bipolar nebular shell visible in the inset of Fig. 3. The rim appears 
fainter in the middle in all lines, except He II A4686 and [Ar IV] 
A4711 - the dust extinction also drops there. The extinction con- 
stant displays some variability with a mean value of 0.183 ± 0.084 
and an rms deviation of 0. 1 1 1 . The [Ar IV] line peaks at the lower 
segment of the rim, whereas He II shows high surface brightness 
spots in three positions along it. Characteristically, these spots are 
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Figure 8. NGC6153 monochromatic line flux (F) and dereddened intensity (I) maps (in units of 10 — 17 ergcm~ 2 s" 1 ), along with the extinction constant 
employed, c(H/3). Blank spaxels: three in 2nd row from bottom (broken Argus fibres); two per corner (fibres reserved for sky subtraction). In this and 
subsequent figures spaxels (3-5, 4) suffer from large uncertainties. The central star position is at coordinates (9, 3). North-west is up (cf. Fig. 2). 



all shifted by one spaxejf] towards the direction of the PN nucleus, 
relative to local enhancements of the H/3 emission, and should cor- 
respond to higher ionization zones perhaps associated with the star- 
facing side of an expanding shell, that is, zones that lie deeper in the 
nebula from the edge of the main ionization front. The C II A4267 
and O II A4649 ORL maps have similar appearances and the lines 
also peak at the lower part of this nebular boundary. 



The seeing was less than 1.5 spaxels wide during these observations. 



4.2 Plasma physical conditions 

The dereddened emission line maps discussed above were used to 
compute line ratio maps and hence the electron temperature (T e ) 
and density (iV e ), and heavy ion abundances in the PNe relative to 
hydrogen. The nebular to auroral [O III] A4959/A4363 ratio was 
used to derive T e , and the [Ar IV] A4711/A4740 ratio was used to 
derive N e . Atomic data for [O III] and [Ar IV] were taken from 
Mendoza (1983) and Zeippen et al. (1987) respectively. Mean val- 
ues of dereddened ratios of lines used in this study, and the cor- 
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Figure 9. NGC7009 monochromatic line flux (J?) and dereddened intensity (/) maps (in units of 10~ 17 erg cm -2 s —1 ), along with the extinction constant 
employed, c(H/3). Blank spaxels: three in 2nd row from bottom (broken Argus fibres); two per corner (fibres reserved for sky subtraction). Spaxels (3-5, 4) 
typically suffer from large uncertainties. The central star position is at coordinates (3, 12). North-east is up (cf. Fig. 3). 



responding plasma conditions for the PNe are reported in Table 2. 
These values are in very good agreement with results from the long- 
slit surveys. 

Table 2 also tabulates the fractional mean-square tempera- 
ture fluctuation parameters, t A , and average electron temperatures, 
To, a, across the IFU area for each PN (that is, in the plane of the 
sky). We have used equations 11 and 12 of Mesa-Delgado et al. 
(2007) for the derivation of these quantities (see also Rubin et al. 
2002), based on the observed (T e , N c ) spatial distributions. Very 
small values of the t 2 A parameter in the plane of the sky are found 



in all three nebulae. The large uncertainties attached to (t A , To,a) 
for NGC6153 are dominated by the relatively large uncertainties 
in the per spaxel determination of 7V C for this nebula. In the case of 
NGC 7009 our t\ estimate of 0.00075 is about five times smaller 
than that derived by Rubin et al. (2002) from //5T/WFPC2 imagery 
of the whole nebula in the [O III] lines. These results do not pre- 
clude the existence of large scale temperature variations along a 
given line of sight. For instance, Barlow et al. (2006) have pre- 
sented preliminary evidence for large line-of-sight T c ([0 III]) vari- 
ations in NGC 6153; however, as discussed in previous works, such 
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Table 2. Mean dereddened line ratios, physical conditions and fractional temperature fluctuation parameters. F(H/3) are observed integrated fluxes across the 
IFUs. The line intensities are in units of /(H/3) = 100. See text for details. 





NGC5882 


NGC6153 


NGC 7009 


c(H/3) 


0.45 ± 0.24" 


1.15 ±0.12 


0.18 ±0.11 


F(H/3) (ergcm- 2 s" 1 ) 


2.66 ■ 10- 11 


3.12 • 10- 12 


3.77 ■ 10- 11 


/(He I A4471) 


5.47 ± 0.23 


6.17 ±0.34 


4.90 ± 0.48 


/(He ii A4686) 


4.67 ± 3.27 


17.4 ± 8.68 


25.60 ± 8.27 


/(C II A4267) 


0.491 ± 0.071 


2.95 ± 0.58 


1.04 ±0.41 


/(O II A4089) 


0.203 ±0.108 


0.679 ±0.173 


0.343 ± 0.098 


/(O II A4649) 


0.433 ± 0.066 


1.88 ±0.41 


0.882 ±0.152 


/([Ne III] A3967) 


30.09 ± 1.77 


32.70 ± 1.68 


44.15 ±6.71 


/([O III] A4959) 


342.3 ±11.2 


291.6 ± 14.2 


372.3 ± 25.4 


/([Ariv] A4711) 


1.97 ± 0.33 


2.66 ± 0.75 


5.04 ± 0.88 


[ArlV] A4711/A4740 


0.867 ± 0.085 


0.972 ± 0.086 


0.896 ± 0.034 


[0 in] A4959/A4363 


250 ±5 


274 ± 47 


187 ± 18 


AT e (cm" 3 ) 


5160 ± 1200 


3650 ± 1400 


4780 ± 490 


To(K) 


9520 ± 190 


9270 ± 410 


10 340 ±380 


T ,A (K) 


9560 ± 2260 


9180 ± 11290 


10270 ± 1970 




0.00034 ± 0.00021 


0.00180 ±0.00234 


0.00075 ± 0.00016 



a These are the uncertainties associated with the variation of the quantity across the IFU region (the rms deviation from the mean). They cannot be directly 
compared to the formal error on an integrated flux measurement for the same region: the uncertainties on the total integrated fluxes and ratios are > 10 times 
smaller. 



large scale temperature gradients cannot be the main cause of the 
abundance discrepancy problem (cf. Section 1). Furthermore, Bar- 
low et al. (2006 and in preparation) discussed new evidence that 
O II ORLs in NGC 6153 and 7009 arise from cool plasma, based 
on their narrower line widths compared to that of the [O III] A4363 
line. 

4.2.1 NGC 5882 

In Fig. 10 maps for NGC 5882 of the dereddened [Ar IV] and 
[O III] diagnostic ratios are presented, along with their respective 
uncertainties determined from the line fitting (which incorporate 
the uncertainties derived from the pipeline processing), and the cor- 
responding nebular N e and T maps and their errors. The uncer- 
tainties on the [Ar IV] ratio are of the order of 10 per cent, whereas 
those on the [O III] ratio are typically smaller. The uncertainties on 
the N c map are of the order of 10 per cent or less along the bright 
PN shell, but increase in lower surface brightness regions (towards 
the centre, and towards the south-east), where the [Ar IV] lines are 
weaker. The mean density is 5160 ± 465 cm -3 with an rms devi- 
ation of 1200 cm- 3 . The T c uncertainties are of the order of five 
per cent or less. The mean temperature is 9520 ± 50 K with an rms 
deviation of 190 K. The eastern part of the PN has the lowest (T c , 
A?c). The observed range in temperature is ~ 1000 K; a couple of 
local maxima are seen, the most conspicuous of which coincides 
with the H/3-bright patch on the northwest part of the nebular shell. 
The mean electron temperature and density in the masked halo re- 
gion, contained in the lower part of the IFU aperture [86 spaxels 
excluding (3^, 4)], are somewhat lower than in the main shell, at 
9380 ± 150 K and 2530 ± 600 crrr 3 respectively. 

4.2.2 NGC 6153 

The N c and T a maps of NGC 6153 (with the accompanying T c un- 
certainties only) are presented in Fig. 11 (top). This PN shows a 
temperature range across the IFU of ~ 1800K, with a mean of 
9270 ± 150 K and an rms deviation of 410 K. The demarcation line 



running diagonally across the error map is caused by the transi- 
tion from single Gaussian profile fits to the [O III] AA4363, 4959 
lines (lower left) to double peaked profiles (upper right). The errors 
are typically larger across the transition zone because of the in- 
creased uncertainty of fitting marginally resolved lines with double 
profiles. Where double Gaussian profiles were fitted the combined 
flux of both was used in the analysis. The spectral resolution was 
high enough to allow a separate T c determination for both nebular 
velocity components but, as the line S/N ratio was not sufficiently 
high, this was not attempted. 

There is a distinct negative temperature gradient from the nu- 
cleus towards the outer nebula (with increasing Y spaxel coor- 
dinates out to an offset of 10 arcsec from the nucleus; see also 
Fig. 15b). The highest values are observed near the nucleus, i.e. 
around spaxel (9, 3), but also in the neighbourhood of (2, 7). A T c 
gradient of very similar magnitude along the minor axis of the PN 
was also observed by Liu et al. (2000), whose fixed long-slit po- 
sition is contained within the IFU orientation (the IFU long-axis 
covers roughly half the extent over which those authors performed 
their spatial analysis). Our observations further resolve distinct lo- 
cal minima such as those at ~ (13, 8). Liu et al. concluded that the 
major part of the T c gradient cannot be due to contamination of the 
A4363 line from recombination excitation. In Sections 5.4 and 6 the 
temperature gradient is discussed in the light of the inhomogeneous 
dust distribution and the higher dust extinction in the inner nebula. 

From our analysis the N e surface variation in NGC 6153 does 
not show a clear correlation with distance from the nucleus. The 
range in [Ar IV] density is ~ 5000 cm -3 with a mean of 3650 ± 
1520 cm -3 and an rms deviation of 1400 cm -3 . There is some 
evidence for slightly higher densities over the H/3 bright patch, 
and lower densities at spaxels beyond its outer edge [at ~ (3, 19); 
bottom-left IFU region). A similar rise in density was also seen by 
Liu et al. roughly correlating with H/3 surface brightness. No firm 
correlation between H/3 and N e can be however established from 
the present analysis as the density uncertainties per spaxel are fairly 
high. 
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Figure 10. NGC5882: maps of the dereddened [Ar IV] and [O III] line ratios and respective errors and of the derived electron density and temperature (and 
respective errors). The density map shows overplotted F(H/3) isocontours. 



4.2.3 NGC 7009 

This nebula shows a mean T c of 10 340 ± 90 K with an rms devia- 
tion of 380 K and a range across the IFU of 1700 K. In contrast to 
NGC 6153, the temperature in NGC 7009 (Fig. 11, bottom) is not 
the highest close to the nucleus at map coordinates (3, 12), but at a 
position ~ 5 arcsec north from it at (12, 13). That position coincides 
with the He II-bright high excitation zone discussed in Section 4.1.3 
(where the He 2+ /He ionization fraction shows a local maximum; 



see also Fig. 14, bottom row), and marks the edge of the ionization 
front of the inner nebular bubble. The N c surface variation across 
the IFU shows enhancements in the top and centre-bottom regions, 
approximately corresponding to local H/3-bright spots; it is lower 
over spaxels (14, 13-20) where the H/3 surface brightness is low. 
The mean electron density across the IFU is 4780 ± 320 cm -3 with 
an rms deviation of 490 cm~ 3 and a range of ~ 3000 cm -3 . 
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Figure 11. Electron density and temperature maps for NGC6153 (top) and NGC7009 (bottom). The density maps show overplotted F(H/3) isocontours. 
Spaxels (3-5, 4) are not considered in the analysis. 



4.3 Chemical abundances 

In this section maps of the chemical abundances for the PNe are dis- 
cussed. The 2+ /H + abundance ratio for all nebulae was derived 
from the [O III] A4959 forbidden line, and also from the O II A4089 
3d-4f and A4649 3s-3p ORLs. The C 2+ /H+ratio was derived from 
the C II A4267 ORL. The He+/H+and He 2+ /H+ abundance ratios 
were derived from the He I A4471 and He II A4686 recombina- 
tion lines respectively; the sum of the two yielded total He/H abun- 
dance ratios for the three nebulae. The T c ([0 III]) and iV ([Ar IV]) 
maps discussed above were adopted in all cases. The reader is re- 
ferred to Tsamis et al. (2003b, 2004) for details on the method of 
calculating ionic abundance ratios based either on heavy element 
CELs or helium and heavier element ORLs. Effective recombina- 
tion coefficients are taken from: Storey (1994) for O II A4649; Liu 
et al. 1995 for O II A4089; Davey, Storey & Kisielius (2000) for 
C II A4267 (including both radiative and dielectronic processes); 
Storey & Hummer (1995) for He I and He II. Table 3 summa- 
rizes the mean values from the maps presented in Figs. 12-14, 
including the abundance discrepancy factors for 2+ . For all tar- 
gets, the mean ionic C 2+ /H+ and 2+ /H+ ORL abundance ratios 
are higher than the corresponding solar abundances of carbon and 
oxygen (e.g. Lodders 2003), by factors of ~ 10 (NGC6153), ~ 
3 (NGC7009), and ~ 2 (NGC5882). We posit that these over- 
abundance factors point towards the presence of heavy element- 
rich (hydrogen-poor) regions within the nebulae. The total carbon 



and oxygen recombination-line abundances, relative to hydrogen, 
for these regions should be even higher. The He/H abundances are 
also higher than solar for all PNe (by up to 40 per cent in the case of 
NGC6153). In the presented maps the posited hydrogen-deficient 
zones are identified as undulations in the ORL-based diagnostics. 



4.3.1 NGC5882 

Abundance maps for NGC 5882 are presented in Fig. 12; the cor- 
responding error maps are shown for the oxygen diagnostics. The 
2+ /H + ratio based on the [O III] A4959 CEL has a mean value 
of (4.28 ± 0.02) ■ 10 -4 with an rms deviation of 2.97 ■ 10" 5 . 
The abundance of this ion shows a local maximum at the position 
of the PN nucleus at map coordinates (6, 15), surrounded by a ~ 
19 arcsec 2 area where the abundance is depressed, and two further 
maxima at positions symmetrically offset by ~ 2 arcsec from the 
nucleus on the north-south direction. The CEL abundance of this 
ion is further depressed over the bright nebular shell but increases 
over the mean value on the eastern side of the PN, where the H/3 
surface brightness and T c ([0 III]) are low. In the masked halo re- 
gion the 2+ /H+ CEL ratio is (4.52 ± 0.28) ■ 10" 4 . 

The 2+ /H + abundance ratio based on the O II A4649 3s-3p 
ORL shows a mean of (9. 10 ± 1 .60) ■ 10~ 4 with an rms variation of 
1.22 ■ 10 -4 . The lowest values of the ratio are found in the south- 
east corner of the IFU. The ionic abundance based on this line is 
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Table 3. Mean chemical abundances and oxygen abundance discrepancy factors (the ratio of 2 + ORL over CEL determinations). The listed uncertainties are 
the formal errors on the mean. The rms deviations from the mean (associated with the variation of the quantity across the IFU region) are quoted in the text. 







NGC 5882 


NGC 6153 


NGC 7009 


2 +/H+ 


[O III] A4959 


(4.28 ± 0.02) ■ 10- 4 


(4.09 ± 0.26) • 10- 4 


(3.54 ± 0.03) ■ 10- 4 


2 +/H+ 


O II A4089 


(1.83 ±0.21)- 10~ 3 


(6.12 ± 0.63) • 10~ 3 


(3.12 ±0.25)- 10~ 3 


ADF(0 2 +) 


A4089/A4959 


4.33 ± 0.49 


16.1 ±2.0 


9.10 ±0.74 


2 +/H+ 


O II A4649 


(9.10 ±1.60) ■ 10- 4 


(4.17 ±0.61)- 10- 3 


(1.78 ±0.13)- 10" 3 


ADF(0 2 +) 


A4649/A4959 


2.14 ±0.35 


10.9 ±2.3 


5.14 ±0.35 


C 2 +/H+ 


C II A4267 


(4.68 ± 0.37) • 10~ 4 


(2.81 ± 0.17) ■ 10~ 3 


(1.01 ±0.05) • 10~ 3 


He+/H 


He I A4471 


0.1042 ±0.0023 


0.1182 ±0.0062 


0.0950 ± 0.0040 


He 2 +/H 


He II A4686 


0.0039 ± 0.0001 


0.0142 ± 0.0010 


0.0213 ± 0.0003 


He/H 


(He+ +He 2 +)/H 


0.1081 ±0.0023 


0.1324 ±0.0063 


0.1163 ± 0.0040 



perhaps slightly depressed over a four spaxel area in the vicinity 
of the nucleus with a value of (6.00 ± 3.13) • 10~ 4 . On the other 
hand, the 2+ /H+ ratio based on the O II A4089 3d-4f line shows 
a very different surface distribution relative to the A4649 transi- 
tion. Whereas the mean value is (1.83 ± 0.21) ■ 10~ 3 with an rms 
deviation of 9.76 • 10~ 4 , the region in the vicinity of the PN nu- 
cleus over a nine spaxel area shows a three times higher abundance 
with a value of (5.52 ± 0.61) • 10~ 3 . From the long slit surveys 
it was established that, on average, the A4089 transition typically 
returns higher O + /H + ratios than the A4649 transition when both 
line emissivities are calculated under the assumption of T e ([0 III]) 
for their emitting regions. This effect, when attributed to the differ- 
ent T e dependence of the emissivities of the two classes of O II tran- 
sitions (3s-3p versus 3d-4f), has been used to provide evidence for 
the existence of cold plasma regions in several PNe consistent with 
the existence of hydrogen-deficient regions (Tsamis et al. 2004; Liu 
et al. 2004). In Tsamis et al. (in preparation) we will revisit this is- 
sue and present a detailed analysis of the O II A4089/A4649 ratios 
from the current data-set based on the latest effective recombina- 
tion coefficients. In the masked halo region the 2+ /H+ ORL ratio 
from the A4649 and A4089 lines is (6.96 ± 1.21) • 10~ 4 and (1.37 
± 0.49) • 10~ 3 respectively, rather consistent with the values across 
the main shell of the nebula. 

We also present the respective 2+ ADF maps, that is 
the 2+ A4649/0 2+ A4959 and 2+ A4089/O 2+ A4959 quantities. 
These appear different (caused by the different behaviour of the 
2+ ORL maps discussed above) with the former showing local 
maxima along the nebular shell and the latter showing a sharp max- 
imum near the PN nucleus. 

The C 2+ /H+ ratio derived from the C II A4267 ORL shows 
a mean of (4.68 ± 0.37) • 10~ 4 with an rms deviation of 6.82 
• 10 -5 . Higher values are observed at the northwest and south- 
east segments of the H/3-bright shell. The value of the ratio over 
four spaxels in the vicinity of the central star is (4.26 ± 1.11) • 
10~ 4 and is therefore consistent with the mean across the IFU (the 
strong stellar continuum in the neighbourhood of the PN nucleus 
can sometimes introduce higher errors in the line fitting - higher 
spectral resolution data would help in such cases). In the masked 
halo region this abundance ratio is (3.88 • 0.42) • 10~ 4 . 

Finally, Fig. 12 (bottom) shows the helium abundance maps 
of NGC 5882. The mean He/H ratio is 0.1081 ± 0.0023 with an 
rms deviation of 0.0029. The He + /H + abundance ratio is depressed 
in the central and western nebular regions, but is enhanced in the 
south and eastern area which is bordering on the PN halo. In con- 
trast, the He 2+ /H + abundance ratio shows two prominent maxima 
offset symmetrically on either side of the PN nucleus along the 



east-west axis. The He /H ratio is inversely correlated with the 
He 2+ /H+ratio across the IFU region (see Table Al). In the masked 
halo region the He+/H+ and He 2+ /H+ ratios are 0.1 142 ± 0.0056 
and (9.14 ± 6.10) • 10~ 4 respectively. 



4.3.2 NGC 6153 

Abundance maps for NGC 6153 are presented in Fig. 13. The 
2+ /H+ ratio based on the [O III] A4959 CEL has a mean value 
of (4.09 ± 0.26) ■ 10~ 4 with an rms deviation of 8.51 • 10~ 5 . The 
forbidden line abundance of this ion shows an inverse correlation 
with T ([0 III]), being high where the plasma temperature is low 
(characteristically in the outer nebula mainly; see Section 5) and 
vice-versa. In contrast, the 2+ /H + ratio based on the O II AA4089, 
4649 ORLs generally increases towards the inner nebula, being 
proportional to the 1(0 Il)/J(H/3) brightness ratio surface distri- 
bution. Both of the O II recombination line abundance diagnostics 
show local maxima in the vicinity of the PN nucleus, but as in the 
case of NGC 5882 (and for the same reasons), their respective maps 
have a different appearance. For this reason, the A4089/A4959 ADF 
map shows values as high as ~ 40 and the A4649/A4959 ADF map 
values as high as ~ 20 at spaxels near the PN nucleus. 

In the same vein, the C 2+ /H + ratio derived from the C II 
A4267 ORL, with a mean of (2.81 ± 0.17) • 10~ 3 and an rms de- 
viation of 5.72 • 10~ 4 , strongly peaks near the PN nucleus as well. 
Characteristically, both the ionic carbon abundance and the ionic 
oxygen abundance (from the O II A4649 line) peak in the neigh- 
bourhood of spaxel (10, 2) about 2 arcsec from the nucleus; the 
2+ /H+ ratio from A4649 has a value of (6.14 ± 0.81) ■ 10~ 3 over 
a three spaxel area when the mean over the IFU is (4.17 ± 0.61) • 
10~ 3 ; this represents a local enhancement of 47 per cent over the 
mean. Across the same region the C 2+ /H + ORL abundance is (3.92 
± 0.27) ■ 10~ 3 , that is 40 per cent higher than the mean. It there- 
fore seems that in the case of NGC 6153 the PN nucleus, either in 
itself or via its hard radiation field, is strongly implicated in the 
causes of the abundance discrepancy, and in the high abundances 
derived from oxygen and carbon ORLs. The peaking of the C II 
and O II abundance diagnostics in the central PN regions had also 
been inferred from the long-slit analysis (Liu et al. 2000). 

In Fig. 13 (bottom) we show the helium abundance maps of 
NGC 6153. The mean He/H ratio is 0.1324 ± 0.0063 with an rms 
deviation of 0.0039. The He + /H + abundance ratio is lower in the 
inner nebular regions (mid-top IFU area), but is enhanced over its 
mean value in the outer nebula. In contrast, the He 2+ /H + ratio is 
higher in the inner PN (see also Section 5). As with NGC 5882, 
there is an inverse linear correlation between the He + /H + and 
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Figure 12. NGC5882: 2 +/H+ abundance ratios from [O III] forbidden and O II recombination lines, C 2 +/H+ abundances from the C II A4267 ORL [with 
-F(H/3) isocontours overplotted], the corresponding oxygen abundance discrepancy factors, and ionic and total He/H abundances from the helium ORLs. 



He 2+ /H + ratios across all spaxels (see Table Al). The total He/H 
abundance possibly shows a small positive gradient towards the in- 
ner nebula (top-left IFU area). 

4.3.3 NGC 7009 

In Fig. 14 abundance maps for NGC 7009 are presented. The 
2+ /H+ ratio based on the [O III] A4959 CEL has a mean value 
of (3.54 ± 0.03) • 10~ 4 with an rms deviation of 4.42 • 10" 5 . In 



this PN also, the forbidden-line abundance of the O + ion shows an 
inverse correlation with T c ([0 III]), being high where the plasma 
temperature is low (see Section 5) and vice-versa. The abundance 
is generally depressed over a triangular region in the middle of the 
IFU. In contrast, over the same area the O II and C II recombina- 
tion line abundance diagnostics are enhanced. There are prominent 
maxima in the 2+ /H + ORL ratio in the vicinity of the PN nucleus 
near map coordinates (3, 12) where the ratio based on the A4649 
line has a value of (3.36 ± 0.62) ■ 10 -3 over a two spaxel area; this 
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Figure 13. NGC6153: 2 +/H+ abundance ratios from [O III] forbidden and O II recombination lines, C 2+ /H+ abundances from the C II A4267 ORL 
[with _F(H/3) isocontours overplotted], the corresponding oxygen abundance discrepancy factors, and ionic and total He/H abundances from the helium ORLs. 
Spaxels (3-5, 4) are not considered in the analysis. 



constitutes a local enhancement of 89 per cent over the correspond- 
ing mean across the map of (1.78 ± 0.13) • 10~ 3 . The respective 
ADF maps show peaks near the nucleus as well. The C 2+ /H + ra- 
tio peaks close to the nucleus with a value of (1.14 ± 0.16) ■ 10 -3 
over a four spaxel area, but also near the edge of the H/3-bright rim, 
some 5 arcsec away from the central star; it has a value of (1.49 ± 
0.04) • 10 -3 at coordinates (11, 18). Its mean value across the IFU 
is (1.01 ± 0.05) ■ 10~ 3 with anrms deviation of 1.57 ■ 10" 4 . 



In Fig. 14 (bottom) the helium abundance maps of NGC 7009 
are shown. The mean He/H ratio is 0.1324 ± 0.0063 with an rms 
deviation of 0.0039. The He + /H + abundance ratio is generally 
low over the triangular area where the 2+ and C 2+ ORL abun- 
dances are enhanced, whereas the He 2+ /H + ratio is correspond- 
ingly high there, peaking near the rim of the inner nebular bubble. 
There is a strong linear inverse correlation between the He + /H + and 
He 2+ /H + ratios in this nebula too (see Table Al). The total He/H 
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Figure 14. NGC7009: 2 +/H+ abundance ratios from [O III] forbidden and O II recombination lines, C 2+ /H+ abundances from the C II A4267 ORL 
[with _F(H/3) isocontours overplotted], the corresponding oxygen abundance discrepancy factors, and ionic and total He/H abundances from the helium ORLs. 
Spaxels (3-5, 4) are not considered in the analysis. 



abundance shows slightly lower values in the top of the IFU map, 
but does not vary across the edge of the inner nebular bubble where 
the He 2+ /He ionization fraction peaks. 



5 CORRELATIONS 

In this section we investigate several correlations between physical 
properties of the PNe that can be established following this analy- 
sis. Evidence for a few of these existed from the long-slit surveys, 
whereas others are new. A unique feature of this study, however, is 
that the large number of data points (~ 300) per PN allowed us to 
put all these on a firm footing for the first time. Moreover, the 2D 
aspect of this analysis is particularly revealing, and in the follow- 
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ing paragraphs it shall be brought to the fore. We present results for 
each nebula individually so as to first highlight the differences and 
then address the similarities; NGC 6153 is examined first. 

5.1 NGC 6153 

The ADF and the [O III] electron temperature for this PN show 
a very distinct correlation with distance from the PN nucleus. We 
have taken the Argus maps of these physical quantities (Figures 1 1 
and 13) and binned all non-zero value spaxels at equal radii out- 
wards from the position of the central star (with a radial increment 
of one spaxel); in Fig. 15a-b the result is shown. As the [O III] tem- 
perature rises from the outer to the inner regions of the PN from ~ 
8400 to 10 000 K, characteristically the oxygen ADF rises as well. 
The quantity plotted in Fig. 15a is the A4649/A4959 ADF (and is 
also used in Figs. 16-18), but a linear correlation is obtained for 
the A4089/A4959 ADF as well. 

Fig. 15c shows the relationship between the oxygen ADF and 
T c ([0 III]) for the 296 spaxels with a measurement of these quanti- 
ties. The solid line is an error-weighted least-squares fit to the data 
with a correlation coefficient (r) of 0.92. In Fig. 15d the inverse 
relationship between the 2+ /H + abundance ratio based on the 
[O III] A4959 line and T c ([0 III]) is shown; the correlation is very 
tight in this case too with a linear correlation coefficient of —0.98. 
In Figs. 15 and 16, black solid lines are least-squares fits weighted 
by the errors, red circles denote data from spaxels with coordinates 
(from the corresponding Argus maps) of (X, Y) = (1-14, 1-11), 
whereas blue diamonds correspond to coordinates (1-14, 12-22). 
This colour scheme allows us to examine the behaviour of nebular 
properties in nebular regions that are near and far out from the cen- 
tral star (red and blue symbols respectively; cf. also the orientation 
of the Argus IFU long axis viz. the position of the PN nucleus in 
Fig. 2). From the juxtaposition of Fig. 15a-b and 15c-d a consis- 
tent picture emerges. Regions closer to the nucleus of NGC 6153 
have higher temperature, lower forbidden line 2+ /H + abundance 
and exhibit higher 2+ ADFs than regions farther out in the neb- 
ula. Some unavoidable intermixing of the colour-coded data is seen 
due to projection (line of sight) effects and/or due to the fact that 
the steady progression from lower to higher temperatures has only 
been highlighted with a two-colour scheme. 

In Fig. 15e-f the correlations between the ORL abundance 
ratios C 2+ /H+ and 2+ H + versus T e ([0 III]) are plotted; the cor- 
relation coefficients of the error-weighted fits are in this case 0.88 
and 0.70 respectively. These plots show that the recombination 
line abundances of doubly ionized carbon and oxygen, which are 
temperature-insensitive to a large degree, (i) rise towards the cen- 
tral nebular regions, and (ii) are typically the highest where the for- 
bidden line [O III] temperature is highest. The 2+ /H + ORL ra- 
tio used in the correlations discussed and plotted in this section is 
based on the O II A4649 line. 

We also find a correlation of the C 2+ /H + ORL abundance ra- 
tio and the He 2+ /He ionization fraction versus the extinction con- 
stant, c(H/3); we performed least-squares fits for the inner nebula 
(~ 148 data points) and show the results in Fig. 15g-h. Fig. 15i 
shows that an inverse trend is established between the oxygen ADF 
and the dereddened H/3 flux from each Argus spaxel. The trend is 
clearly non-linear. Regions of low surface brightness (in H/3) that 
lie in the inner nebula exhibit the highest oxygen ADFs, whereas 
the low ADF 'tail' corresponds to brighter zones farther out in the 
nebula (most points in this tail correspond to spaxels of the last 
few IFU rows - large Y coordinates - on the far side from the PN 
nucleus). 



The correlations between the ORL abundance ratios C /H 
and 2+ /H + versus the forbidden-line 2+ /H + abundance ra- 
tio are shown in Fig. 16a-b, where an inverse trend is seen 
in both cases with r's of —0.86 and —0.69 respectively. These 
are mainly driven by the very strong linear correlation be- 
tween the temperature-dependent forbidden-line 2+ /H + ratios 
and T c ([0 III]) that was discussed earlier. In stark contrast, a posi- 
tive correlation between the temperature-insensitive ORL C 2+ /H + 
and ORL 2+ /H + abundance ratios is established with r — 0.87 
(Fig. 16c). The fact that the trends between the ORL-CEL and 
ORL-ORL abundance diagnostics lie in completely opposite direc- 
tions shows that (i) vastly different assumptions govern the physical 
origins and the metallicity of the gas traced by ORL diagnostics 
on the one hand and CEL diagnostics on the other, and that (ii) 
consistent nebular metallicity estimates can in principle only be 
derived using pure ORL/ORL or CEUCEL emission line ratios^ 
This further tells us that the C 2+ /0 2+ (pure ORL) abundance ra- 
tio is approximately constant throughout the portion of NGC 6153 
observed in this programme with a mean value of 0.68 ±0.11. 

Fig. 16d-e shows that the C 2+ /H + abundance ratio and the 
oxygen ADF are positively correlated with the He 2+ /He ioniza- 
tion fraction with r's of 0.85 and 0.86 respectively. The He 2+ /He 
fraction rises close to the PN nucleus where the ionization degree 
of the nebular plasma is higher due to the harder radiation field, so 
these plots are again consistent with an origin for the ORL-emitting 
heavy element-rich plasma in the inner regions of the nebula. On 
the other hand an inverse linear trend is established between the 
2+ /H+CEL abundance ratio and He 2+ /He (r = -0.87; Fig. 16f). 

In Fig. 16g-h the correlations between the C 2+ /H + ORL 
abundance, the oxygen ADF, and the 2+ /H+ ORL abundance 
versus the He/H abundance ratio are shown. Positive, but weaker, 
trends are seen in all cases with r's of 0.68, 0.59, and 0.55 respec- 
tively. These plots also reveal that typically regions in the inner neb- 
ula show a somewhat higher helium abundance (but note that this 
large He/H is challenged in Section 6). Finally, we do not find any 
obvious trends of the C 2+ /H + or 2+ /H + ORL abundances (or of 
the 2+ ADF) with electron density. The detailed parameters of our 
regression analysis for this nebula are listed in Table Al. Parame- 
ters of fits obtained when the 2+ /H + ORL abundance is based on 
the A4089 line are also listed in Table Al (the corresponding plots 
are not shown): the fits are all in the same direction as previously, 
and of similar strength, but with a slightly larger scatter [as the per 
spaxel uncertainties on J(A4089) are somewhat higher]. 

5.2 NGC 7009 

In this section we examine whether the various correlations be- 
tween diagnostics that were established for NGC 6153 hold for this 
nebula too. In order to check for any markedly different behaviour 
between nebular regions the data from the Argus maps have been 
separated according to the following scheme: in the following plots, 
data from spaxels with coordinates - see Fig. 9 - (X, Y) = (1-9, 5- 
19) are represented by red circles, whereas data points outside this 
rectangular area are marked by blue diamonds. The red symbols 
thus correspond to a nebular area of lower mean surface bright- 
ness (in dereddened H/3) and of somewhat closer proximity to the 
PN nucleus than the blue symbols (cf. Fig. 3). The blue symbols 
comprise the zones of higher surface brightness and specifically 
the bright H/3 rim of the inner PN shell which runs from top to 

5 This may not work for helium however; see Section 6. 
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Figure 15. NGC6153: Correlations between various physical quantities: (a) The 2 + ADF and (b) the [O III] T c radial profiles resulting from binning the 
respective Argus maps at equal radii from the central star. The full radial offset from the PN nucleus (20 spaxels) is 10.4 arcsec; (c)-(i) Red circles correspond 
to Argus IFU spaxel coordinates (X, Y) = (1-14, 1-1 1) and blue diamonds to coordinates (1-14, 12-22) - see Fig. 8. The solid lines are weighted fits; (g)-(h) 
The solid lines are weighted fits to the red data points only. See text for details. 



bottom in the H/3 image of the nebula (cf. Fig. 9 and the inset of 
Fig. 3 where this sharp boundary is clearly seen in the HST im- 
age). This nebula is however larger than NGC6153 and the Argus 
IFU targeted the central regions of the inner bipolar shell; hence 
the adopted grouping of spaxels in the plots that follow may not 
be as clear-cut regarding their distance from the nucleus as it was 
for NGC 6153 where the IFU was both more favourably placed and 
captured a larger portion of the PN. 

In Fig. 17a the inverse relationship between the 2+ /H + 
abundance ratio based on the [O III] A4959 line and T e ([0 III]) 
is plotted for the 296 spaxels with a measurement of these quanti- 
ties. The solid line is an error-weighted least-squares fit to the data 
with a correlation coefficient of —0.94. The positive trend between 
the oxygen ADF and T c ([0 III]) is shown in Fig. 17b and has a 
correlation coefficient of 0.78. 

In Fig. 17c-d the positive correlations between the ORL abun- 
dance ratios C 2+ /H+ and 2+ /H + versus T c ([0 III]) are shown; 
the r's of the error-weighted fits through the 135 red data points 
of the inner nebular field are in this case 0.66 and 0.70 respec- 
tively; they are 0.50 and 0.35, respectively, when fitting all the data. 
These plots show that the recombination line abundances of doubly 



ionized carbon and oxygen for NGC 7009 follow the same trend 
with forbidden line temperature as was the case of NGC 6153. The 
2+ /H + ORL ratio used in the correlations plotted in this section 
is based on the O II A4649 line. 

In Fig. 17f-g we plot the correlations between the ORL 
abundance ratios 2+ /H + and C 2+ H + versus the forbidden-line 
2+ H + abundance ratio, where an inverse trend is seen in both 
cases with r's of —0.68 and —0.63 respectively for fits to the red 
data points. When fits to all data are performed the coefficients 
are less than 0.50 due to the larger scatter of some of the blue 
data points which correspond to nebular zones of higher mean sur- 
face brightness. The general behaviour of the trends is the same 
as for NGC 6153. A positive correlation between the temperature- 
insensitive ORL C 2+ /H + and ORL 2+ /H + abundance ratios is 
established (r = 0.86; Fig. 17e). The C 2+ /0 2+ (pure ORL) abun- 
dance ratio is approximately constant throughout the field but there 
are a few (red) points a factor of ~1.27 higher (but within 2a) than 
the mean ratio of 0.57 ± 0.05. 

Fig. 17i-j shows that the C 2+ /H + ORL abundance ratio and 
the 2+ ADF are positively correlated with the He 2+ /He ioniza- 
tion fraction with r's of 0.35 and 0.55 respectively (for fits through 
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Figure 16. NGC 6153: Correlations between various physical quantities with symbols as in Fig. 15; the solid lines are weighted fits. See text for details. 



the red data points only). The 2+ /H + CEL ratios are inversely cor- 
related with the helium ionization state (the weighted fit to all data 
has an r of —0.84). There is more scatter in these plots than was the 
case for NGC 6153: the cluster of low ionization points within the 
dashed squares in Fig. 17i-j (at ADF ~ 5) come from the lower- 
right corner of the corresponding He 2+ /He map and also from IFU 
row Y = 14; they are also of relatively high 2+ /H + (forbidden- 
line) abundance. The main trends identified in these plots are in 
the same vein as for NGC 6153, with high ionization gas exhibiting 
higher oxygen ADFs and higher C 2+ /H + ORL abundance ratios 
than lower ionization, high 2+ /H + CEL regions. 

Fig. 17k-m shows the ORL 2+ /H+ and C 2+ H + abundance 
ratios and the 2+ ADF versus the He/H abundance ratio. Weak 
positive trends are seen in all cases - for the 2+ ADF vs. He/H lin- 
ear fit, r — 0.58. The trends are in the same direction as in the case 
of NGC 6153. We further plot the 2+ ADF versus the dereddened 
H/9 intensity per spaxel in Fig. 17n. No clear trend is seen in con- 
trast to the NGC 6153 case. We do not find any obvious trends of the 
C 2+ /H+ or 2+ /H+ ORL abundances (or of the 2+ ADF) with 
electron density. No firm trend of the heavy element ORL abun- 
dances with the extinction constant is found in contrast to the case 
of NGC 6153. The detailed parameters of our regression analysis 
for this nebula are listed in Table Al. Parameters of fits obtained 
when the 2+ /H + ORL abundance is based on the A4089 line are 



also listed in Table Al (the corresponding plots are not shown): the 
fits are in the same direction as previously. 



5.3 NGC 5882 

This PN is the smallest of the three and its symmetrical inner shell 
was almost fully captured by the Argus IFU. It also displays lower 
mean oxygen ADFs and was meant to be our 'control' target. In this 
section we explore correlations between the same set of quantities 
for NGC 5882 as for the other two nebulae. The Argus spaxels were 
grouped based on their distance from the nucleus in the following 
manner (cf. Fig. 7): points with spaxel coordinates (X, Y) = (4- 
10,11-18) are coded as red circles whereas those exterior to this 
boxed region are coded as blue diamonds. 

Fig. 18a shows the 2+ /H + abundance ratio based on the 
[O III] A4959 line versus T e ([0 III]) for the 203 spaxels with a mea- 
surement of these quantities (that is, excluding the masked halo re- 
gion of the nebula). The solid line is an error-weighted least-squares 
fit to the data with a correlation coefficient of —0.89. The positive 
trend between the oxygen ADF and T e ([0 III]) (Fig. 18b; r = 0.52) 
is not as strong as for the other two PNe. In Fig. 18c-d we plot 
the C 2+ /H + and 2+ H+ ORL abundance ratios versus T c ([0 III]); 
both show loose positive trends (with equal r = 0.27). The 2+ /H + 




Figure 17. NGC7009: Correlations between various physical quantities. The solid lines are weighted fits. Red circles correspond to Argus IFU spaxel 
coordinates (X, Y) = (1—9, 5-19) and blue diamonds correspond to spaxels outside this area (see Fig. 9); (f, g) and (i, j): the solid lines are fits to the red data 
points only. See text for details. 
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Figure 18. NGC 5882: Correlations between various physical quantities. Red circles correspond to Argus IFU spaxel coordinates (X, Y) = (4-10, 1 1-18) and 
blue diamonds correspond to spaxels outside this area (see Fig. 7). The solid lines are linear fits weighted by the errors (only through the blue data points in 
panel 'h'). See text for details. 
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ORL ratio used in the correlations plotted in this section is based 
on the O II A4649 line. 

The 2+ /H+ and C 2+ H+ ORL abundances are plotted against 
the forbidden-line 2+ /H + abundance ratio in Fig. 18f-g; loose 
negative trends are established in both cases, but in the same vein 
as for NGC6153 and NGC7009 where the trends were stronger. 
The C 2+ /H + ORL abundance is positively correlated with the 
2+ /H + ORL abundance ratio (r = 0.61; Fig. 18e) following the 
same trend as for the other targets. The C 2+ /0 2+ (pure ORL) abun- 
dance ratio of NGC 5882 is approximately constant across the IFU 
region with a mean value of 0.52 ± 0. 1 1 . 

A loose negative trend is seen between the 2+ /H + CEL 
abundance and the He 2+ ionization fraction (mostly for low ion- 
ization gas represented by the blue data points with r = —0.68; 
Fig. 18h). The C 2+ /H + ORL abundance shows a flat trend with in- 
creasing He 2+ /He (Fig. 18i), in contrast to the cases of NGC 6153 
and 7009, while the 2+ ADF shows a week positive correlation 
with He 2+ /He (Fig. 18j; r = 0.42). The excitation state of the gas 
in this PN is however the lowest of the three, with less than 10 per 
cent of helium being doubly ionized (compared with values of up 
to 20 and 30 per cent for NGC 6153 and NGC 7009 respectively). 

No clear trends are seen between the doubly ionized carbon 
and oxygen ORL abundances, and the 2+ ADF versus the total 
He/H abundance (Fig. 18k-m) when positive trends were estab- 
lished for the other two nebulae. Finally, no trend is established 
between the 2+ ADF and the H/3 intensity per spaxel either 
(Fig. 18n). The trend between the C 2+ /H+ ORL abundance ratio 
and the extinction constant is flat and is not shown (parameters of 
the fit are however listed in Table Al). Detailed parameters of our 
regression analysis of NGC 5882 are listed in Table Al, including 
those obtained when the 2+ /H + ORL abundance is based on the 
A4089 line; those fits are in the same vein as previously. 

5.4 Correlations: what do they tell us? 

A strong inverse correlation between the CEL-based 2+ /H + 
abundances and T c ([0 III]) is established for all three sources. 
The 2+ /H + CEL abundance decreases by factors of approxi- 
mately 2.0, 1.5 and 1.3 over AT c 's of ~ +1600, +1100 and 
+700 K in NGC 6153, 7009, and 5882 respectively (Figs. 15d, 
17a, 18a). The exponential dependence of CEL-based abundances 
on the adopted electron temperature (the Boltzmann factor) math- 
ematically suppresses most of the scatter in 2+ /H + induced by 
errors in the T c ([0 III]) determination, thus strengthening the cor- 
relation. However, the correlation mostly reflects a physical process 
since the variation of temperature with position is not random (see 
NGC 6153, but also NGC 7009 where regions closer to the nucleus 
have a higher mean temperature), and T e ([0 III]) shows a high de- 
gree of correlation with other essentially independent quantities 
such as the C 2+ /H + and 2+ /H + ORL abundances. A physical 
correlation is expected as well, since [O III] is a dominant coolant 
of the gas: where the concentration of 2+ is lower, T e ([0 III]) 
will tend to be higher. This may be due to a real variation of the 
O/H abundance or to a variation of the 2+ /0 ionization fraction. 
Assuming a constant O/H, the higher forbidden-line temperature in 
the inner nebular regions can be qualitatively understood if part of 
the 2+ ionic component shifts into O s+ ; this is a likely possibil- 
ity since 3+ must coexist with He 2+ (based on ionization stratifi- 
cation considerations), and the 2+ /H + CEL abundance inversely 
correlates with He 2+ /He (Figs. 16f, 17h, 18h). 

Along the same lines, contrasting NGC 6153 with NGC 5882, 
the observed range of T e ([0 III]) across the respective IFU regions 



is larger where the He zone is more prominent (i.e., in the for- 
mer rather than in the latter nebula). The large opacity of He + 
and the relatively large average energy of the corresponding photo- 
electrons both contribute to further heat the high-ionization region. 
That the inner nebula electron temperature can be explained in this 
way should await quantitative analysis by means of photoioniza- 
tion models tailored to these observations. Extra heating can be 
provided by the photoelectric effect on small dust grains that may 
populate the vicinity of the central stars (see Section 6). The variety 
of potential heat sources and the associated uncertainties can make 
potential spatial variations of O/H in the gas phase that emits [O III] 
A4959 difficult to ascertain. One intriguing aspect is the fact that 
T c ([0 III]) continues to decrease outwards in NGC 6153, even in 
regions where He + and O should in principle be minor species. 
It may be that He 2+ can survive relatively farther out from the nu- 
cleus than expected as the gas 'porosity' induced by density fluctu- 
ations can blur out the theoretical ion stratification. 

The second fundamental piece of information from this anal- 
ysis is the clear positive correlation of the C 2+ /H + and 2+ /H + 
ORL abundances versus T c ([0 III]) for all three targets - even for 
the low ADF 'control' NGC 5882. Throughout the IFU region of 
NGC 6153, for instance, these abundances increase by a factor of 
~ 2 (Fig. 15e-f) from the outer to the inner nebula, so that the cor- 
responding increase in the 2+ ADF is a factor of ~ 4 (Fig. 15c). 
This is paradoxical in that ORLs should a priori be enhanced in low 
T e conditions (due to their emissivities obeying an inverse power 
law with electron temperature). If a sizeable fraction of the (optical) 
[O III] flux were to arise from the same parcel of gas emitting a size- 
able fraction of ORL flux, then part of the [O III] emission would 
necessarily be produced at relatively low electron temperatures and 
the (observed) average T c ([0 III]) would be lower (not higher), 
where ORL emission is strongest. Thus, this paradox demonstrates 
beyond reasonable doubt that heavy element CELs and ORLs must 
essentially arise from unrelated gas components. 

The 'dual abundance' picture of two distinct components of 
highly ionized gas at very different temperatures can solve the para- 
dox, but the only known way to maintain a strongly photoionized 
gas at low T is thermostasis by fine-structure (FS) infrared lines. 
This mechanism can only be effective for sufficiently large heavy- 
element abundances (and sufficiently low electron densities so that 
the low critical density FS lines are not quenched). This interpre- 
tation of the data is qualitatively coherent since cool hydrogen- 
deficient ionized gas propitiously emits ORLs. Therefore, at an- 
other quantitative level (separate from ORLs), an essential check 
can be provided by the 2D spatial distribution of infrared FS line 
emission; this should be tractable by means of Spitzer Space Tele- 
scope observations and future Herschel surveys. 



6 DISCUSSION AND CONCLUSIONS 

We have presented the first detailed analysis of integral field unit 
spectrophotometry of three Galactic PNe (NGC 5882, 6153, and 
7009) taken with the 8.2-m VLT and demonstrated the superb ca- 
pability of FLAMES Argus data for the accurate mapping of the 
physical properties of nebulae. Spatially resolved maps of 11.5" 
x 7.2" areas of the PNe with 0.52 2 arcsec 2 spaxels were made 
in the light of [O III] A4959 and [Ne III] A3967 lines, and the H/3 
A4861, He I A4471, He II A4686, C II A4267, O II A4089 and A4649 
recombination lines. The dust extinction and the plasma tempera- 
ture and density were mapped out, along with the ionic abundances 
of helium, C 2+ , and 2+ , relative to hydrogen. The 2+ /H+ ra- 
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tio was derived from the O II optical recombination lines and the 
[O III] A4959 collisionally excited line. Maps of the resulting 2+ 
abundance discrepancy factor were made. The 2+ ADF varies be- 
tween the targets and across the Argus field of view in each case. In 
NGC6153 and NGC7009 the 2+ ADF and the C 2+ ORL abun- 
dance show distinct spikes at or near the nucleus. In NGC 7009 
these quantities further peak at positions coinciding with the high 
excitation boundary between the inner nebular shell and the outer 
envelope; this feature is well registered on both the Argus maps 
and the HST WFPC2 image. In NGC 5882, where the S/N ratio of 
O II ORLs is lower, the 2+ ADF derived from the A4089/A4959 
diagnostics also peaks at the position of the central star. 

In all cases these local maxima are observed over about two 
spaxels (~1 arcsec) and are therefore at the limit of our spatial 
resolution. This means that the hypothetical structures associated 
with the hydrogen-poor plasma should have physical dimensions of 
< 1000 astronomical units (a.u.), assuming a typical PN distance 
of the order of 1 kpc. 

Correlations between the O + ADF and the ionization state of 
the gas as gauged by the He + /He ionization fraction were estab- 
lished for NGC 6153 and 7009; likewise strong correlations of the 
2+ ADF and the C 2+ /H+ and 2+ /H+ ORL abundances versus 
the forbidden line electron temperature were established for these 
nebulae (for NGC 5882 these correlations were weaker but in the 
same direction). This constitutes new evidence that the 'metallic' 
ORLs in these PNe are somehow associated with highly ionized 
plasma, and that heavy element ORLs and CELs arise from largely 
unrelated parcels of gas, in accordance with dual-abundance pho- 
toionization models (Pequignot et al. 2002). We argue that this 
hydrogen-deficient, highly ionized plasma may well (i) have been 
recently ejected from the nucleus in the form of clumped gas, or 

(ii) be part of circumstellar material very close to the star in the 
form of a disk-torus structure (based on the spatial correlation of 
the diagnostics with the position of the central star), or equally that 

(iii) it originates from high excitation gas currently being photoe- 
vaporated from milli-pc scale neutral condensations embedded in 
the nebulae at the limit of our spatial resolution. 

In case (i), the physical origin of the enhanced heavy-element 
ORL emission could be an ensemble of hydrogen-poor clumps anal- 
ogous to the H-poor knots observed in the class of nebulae, such as 
Abell 30 orAbell 58 (e.g. Borkowski et al. 1993; Guerrero & Man- 
chado 1996), which may have undergone a late helium-flash (Iben 
et al. 1983) or could be associated with neon novae (Wesson et al. 
2007). In this context, the range of ADFs found amongst our targets 
could reflect an evolutionary effect with H-poor zones becoming 
progressively more readily observed as the density and tempera- 
ture contrast between them and the ambient expanding nebula in- 
creases over time. One should recall that NGC 5882 is positioned in 
the lower end of the ADF sequence' established from the long-slit 
surveys according to which dense, compact and lower ionization 
PNe exhibit smaller integrated abundance discrepancies than more 
diffuse, larger and higher ionization nebulae. NGC 6153 and 7009, 
on the other hand, appear to be in a more advanced stage of the 
ADF evolution' (Tsamis et al. 2004; Liu et al. 2004). 

In case (ii), the hydrogen-deficient gas could arise from cir- 
cumstellar structures such as those recently discovered by VLT in- 
terferometric observations in association with J-type C-stars (IRAS 
18006-3213 - Deroo et al. 2007a), post-AGB stars (V390 Vel - 
Deroo et al. 2007b) and PNe (Mz 3 - Chesneau et al. 2007). These 
toroidal structures have estimated radii of < 1000 a.u. and there is 
evidence that they contain large amounts of silicate dust. The infer- 
ence about their chemistry is important as the C/O ORL abundance 



ratio derived for our PN sample is lower than unity; this means that 
an association of the hydrogen-poor, O-rich gas with silicate-based 
(i.e. oxygen-rich) dust is not unlikely (see below). 

In case (iii) the H-poor plasma could originate in evaporat- 
ing dusty 'cometary knots ' analogous to those in the Helix nebula 
(NGC 7293; e.g. Meaburn et al. 1992) immersed in high ionization 
nebular zones both near (NGC 5882, 6153, 7009) and farther out 
(NGC 7009) from the stars. The Helix globules however are rich 
in molecular hydrogen (Meixner et al. 2005) and this hypothesis 
would require a mechanism for the efficient removal of hydrogen 
and/or the chemical differentiation of the evaporated gas for the 
spectroscopic signature of H-poor plasma to be observed. Also, in 
the Helix nebula at least, there is no evidence of globules inside the 
high ionization He II zone or of He II emission from the observed 
globules, and these features are mostly associated with emission 
from lower ionization species (O'Dell et al. 2002; O'Dell, Hen- 
ney and Ferland 2007); in contrast, as we have shown, the H-poor 
plasma in our PNe is associated with He II emission. On the other 
hand, this does not preclude that 'molecular' globules can be effec- 
tively destroyed in the inner regions of PNe whilst leaving pockets 
of H-poor plasma behind. This hypothesis can clearly be put to the 
test with observations of the type presented in this work. 

We have presented evidence that, in the central regions of 
NGC 6153, the C 2+ ORL abundance, the He 2+ /He ionization frac- 
tion and the 2+ ADF are weakly correlated with the extinction 
constant. There is marginal evidence that this may hold for the 
inner parts of NGC 7009 as well. This correlation, if real, could 
corroborate either the 'dusty clump evaporation' scenario or the 
'dusty disk' scenario for the origin of the hydrogen-deficient gas. 
Dust-rich regions can contribute to the heating of their plasma envi- 
rons via photoelectric gas heating (e.g. Dopita & Sutherland 2000; 
Stasinska & Szczerba 2001; Weingartner, Draine and Barr 2006) 
and to the metal enrichment of the gas via grain destruction, thus 
bringing about this correlation. According to Dopita & Sutherland 
dust grain photoelectric heating is more effective in regions with a 
population of small grains, it increases with increasing grain metal- 
licity, and is more important in highly ionized regions. This effect 
could also play a role in the observed rise of T c ([0 III]) in the in- 
ner NGC 6153 regions (in parallel to usual photoionization which 
can also enhance the electron temperature in the He 2+ zone), al- 
though coupled dust/plasma models would be needed to ascertain 
and quantify this. Alternatively, since the theoretical H7/H/3 line ra- 
tio decreases with decreasing electron temperature (Storey & Hum- 
mer 1995), and since the extinction constant (based on H7/H/3) 
was computed here adopting T c ([0 III]), the positive correlation 
of c(Hf3) with the ADF may be indicative of the fact that H I lines 
in PNe partly arise from cool gas which emits most of the heavy 
element ORLs, in accordance with dual-abundance models. 

Relatively weak correlations between the 2+ ADF and the 
2+ /H + , C 2+ /H + ORL abundances versus the He/H ratio in 
NGC 6153 and (in the inner regions) of NGC 7009 have been es- 
tablished. This result, which should be backed up by further stud- 
ies, constitutes new observational evidence that corroborates the 
predictions of the dual-abundance PN models of Pequignot et al. 
(2002, 2003). In those models, the heavy element-rich (H-poor) gas 
component, which is responsible for the bulk of the heavy element 
ORL emission, is also helium-rich and as such it contributes non- 
negligibly to the emissivity of the helium ORLs. When this effect 
is not accounted for, as in empirical methods like the one we have 
used in this work, there is an insidious bias towards overestimat- 
ing the helium to hydrogen abundance ratio. Based on our results 
it is therefore clear that total He/H abundances of PNe as tradition- 



24 Y. G. Tsamis et al. 



ally determined from empirical methods, and which assume chem- 
ical uniformity of the plasma, will be systematically overestimated. 
Again detailed modelling in the context of a chemically inhomoge- 
neous nebula is needed to remove this bias. With respect to helium 
it should be noted that, in contrast to He 2+ , He + is not perturbed 
by H-poor material, since nebular regions with high He 2+ /H + give 
the same (He + + He 2+ )/H ratio as regions with low He + /H + . 

In conclusion, our detailed IFS analysis of a few suitable ob- 
jects has allowed us to identify important new elements of the still 
incomplete puzzle of planetary nebula astrophysics, and we have 
put forward testable proposals for the likely nature of some of the 
missing pieces. Furthermore, this study has underlined a number 
of clear requirements that a detailed solution to the 'ORL vs. CEL' 
problem should meet, if it were to be successful. 
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Table Al. Parameters of least-squares regression analysis discussed in Section 5: A - intercept, B - slope, r - correlation coefficient, <r - standard deviation, 
Points - number of data points fitted. Values in parentheses are exponents of base 10. 
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C 2 +/H+ ORL vs. 2 +/H+ A4089 ORL 


3.11(-4) ± 3.31(-5) 


0.230 ±0.011 


0.78 


2.2 


293" 


2 +/H+ A4649 ORL vs. 2 +/H+ cel 


4.78(-3) ± 2.81(-4) 


-8.70 ±0.811 


-0.68 


1.6 


bs 6 


2 +/H+ A4089 ORL vs. 2 +/H+ cel 


6.10(-3) ± 2.53(-4) 


-8.41 ± 0.694 


-0.58 


2.0 


293 a 


C 2 +/H+ ORL vs. 2 +/H+ CEL 


2.15(-3) ± 1.22(-4) 


-3.27 ± 0.350 


-0.63 


1.6 


135 6 


2 +/H+ CEL vs. He 2 +/He 


4.39(-4) ± 2.90(-6) 


-4.22(-4) ± 1.64(-5) 


-0.84 


7.9 


289 


C 2 +/H+ ORL vs. He 2 +/He 


7.39(-4) ± 6.37(-5) 


1.36(-3) ±3.10(-4) 


0.35 


1.9 


135 b 


ADF(0 2 +A4649) vs. He 2 +/He 


1.19 ±0.53 


19.3 ± 2.57 


0.55 


2.6 


135 6 


ADF(O 2 +A4089) vs. He 2 +/He 


5.78 ± 0.28 


15.9 ± 1.62 


0.50 


3.0 


293" 


2 +/H+ A4649 ORL vs. He/H 


-4.79(-4) ± 4.64(-4) 


0.0201 ± 0.0040 


0.28 


2.5 


294 


2 +/H+ A4089 ORL vs. He/H 


-6.38(-3)± l.ll(-3) 


0.0819 ± 0.0096 


0.45 


2.3 


290 a 


C 2 +/H+ ORL vs. He/H 


-3.13(-4)±2.88(-4) 


0.0115 ±0.0025 


0.26 


3.4 


295 


ADF(0 2 +A4649) vs. He/H 


-19.5 ±2.01 


214 ± 17.5 


0.58 


3.3 


295 


ADF(O 2 +A4089) vs. He/H 


-38.5 ± 3.90 


409 ± 34.0 


0.58 


2.9 


293 a 


He+/H+ vs He 2 +/H+ 

lit /ll VS. lit /ll 




_0 700 + 018 

\J.I\jy _L V7.V7 1 


—0.93 


0.86 


293™ 






NGC 5882 








2 +/H+ CEL vs. T c ([0 III]) 


1.71(-3) ±4.52(-5) 


-1.34(-7)±4.72(-9) 


-0.89 


8.2 


204 


ADF(0 2 +A4649) vs. T e ([0 HI]) 


-5.39 ± 0.87 


7.92(-4) ± 9.09(-5) 


0.52 


1.0 


202 


ADF(O 2+ A4089) vs. T e ([0 HI]) 


-21.8 ± 2.62 


2.70(-3) ± 2.74(-4) 


0.58 


1.8 


197 a 


C 2 +/H+ ORL vs. T c ([0 III]) 


-3.55(-4) ± 2.06(-4) 


8.75(-8) ±2.16(-8) 


0.27 


1.8 


203 


2 +/H+ A4649 ORL vs. T ([0 III]) 


-6.16(-4) ± 3.85(-4) 


1.60(-7) ±4.05(-8) 


0.27 


l(-4) 


203 


2 +/H+ A4089 ORL vs. T e ([0 HI]) 


-2.21(-3)± 1.13(-3) 


4.07(-7) ± 1.19(-8) 


0.24 


1.8 


190 a 


C 2 +/H+ ORL vs. 2 +/H+ A4649 ORL 


1.56(-4) ±2.98(-5) 


0.347 ±0.031 


0.61 


1.5 


203 


C 2 +/H+ ORL vs. 2 +/H+ A4089 ORL 


5.13(-4) ± 1.90(-5) 


-0.0164 ± 0.0109 


-0.11 


1.9 


196 a 


2 +/H+ A4649 ORL vs. 2 +/H+ CEL 


1.37(-3) ± 1.12(-4) 


-1.07 ±0.260 


-0.28 


l(-4) 


203 


2 +/H+ A4089 ORL vs. 2 +/H+ cel 


2.95(-3) ± 3.15(-4) 


-3.01 ± 0.739 


-0.28 


1.82 


193 a 


C 2 +/H+ ORL vs. 2 +/H+ CEL 


7.13(-4) ±5.95(-5) 


-0.547 ±0.140 


-0.26 


1.8 


203 
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Table A2. -continued 



Diagnostics 


A 




B 




r 


a 


Points 


2+ /H+ CEL vs. He 2 +/He 


4.47(-4) ± 2.97(- 


-6) 


-8.48(-4) ± 7.62(- 


-5) 


-0.68 


14 


148 c 


C 2 +/H+ ORL vs. He 2 +/He 


4.88(-4) ± 7.70(- 


-6) 


-1.44(-4) ± 1.73(- 


-4) 


-0.06 


1.9 


199 


ADF(0 2 +A4649) vs. He 2 +/He 


1.95 ±0.04 




5.62 ± 0.87 




0.42 


0.31 


200 


ADF(O 2 +A4089) vs. He 2 +/He 


3.52 ±0.13 




12.6 ±3.19 




0.27 


2.7 


200" 


2 +/H+ A4649 ORL vs. He/H 


2.23(-3) ± 2.94( 


-4) 


-0.012 ± 0.003 




-0.30 


l(-4) 


203 


2 +/H+ A4089 ORL vs. He/H 


4.87(-3) ± 1.22(- 


-3) 


-0.029 ±0.011 




-0.18 


2.6 


203" 


C 2 +/H+ ORL vs. He/H 


9.29(-4) ± 1.72(- 


-4) 


-4.12(-3)± 1.59(- 


-3) 


-0.18 


1.9 


203 


C 2 +/H+ vs. c(H/3) 


4.76(-4)±8.1(- 


-6) 


1.38(-5) ± 1.61(- 


5) 


0.06 


1.9 


203" 


ADF(0 2 +A4649) vs. He/H 


7.07 ± 0.84 




-45.4 ± 7.80 




-0.38 


0.31 


203 


ADF(O 2 +A4089) vs. He/H 


12.9 ± 2.46 




-83.1 ±22.7 




-0.25 


2.2 


197" 


He+/H+ vs. He 2 +/H+ 


0.108 ±3.27(-4) 


-1.13 ±0.071 




-0.75 


1.4 


201" 



a Corresponding plot not shown. 

6 Fit only for the inner nebula (coloured red) data points; see text for details. 
c Fit only for the outer nebula (coloured blue) data points; see text for details. 



